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O 

Cn ' Abstract. Using the spectral theory of unitary operators and the the- 

r*^ I ory of orthogonal polynomials on the unit circle, we propose a simple 

C^ . matrix model for the following circular analogue of the Jacobi ensemble: 

n 
l<k<l<n j = l 

with D\e 5 > —1/2. If e is a cyclic vector for a unitary n x n matrix 

m^ • U, the spectral measure of the pair {U, e) is well parameterized by its 

^^ ' Verblunsky coefficients (ao, . . . , a„-i). We introduce here a deformation 

(70, . . . , 7n-i) of these coefficients so that the associated Hessenberg ma- 

rG • trix (called GGT) can be decomposed into a product r-(7o) ■ ■ • r(7„_i) of 

jrt I elementary reflections parameterized by these coefficients. If 70, . . . , 7n-i 

are independent random variables with some remarkable distributions, 
then the eigenvalues of the GGT matrix follow the circular Jacobi dis- 
tribution above. 

fT^ ■ These deformed Verblunsky coefficients also allow to prove that, in 

^ ' the regime 5 = 5{n) with 5{n)/n — >■ /?d/2, the spectral measure and the 

Cn , empirical spectral distribution weakly converge to an explicit nontrivial 

probability measure supported by an arc of the unit circle. We also 

^/ ' prove the large deviations for the empirical spectral distribution. 
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1. Introduction 

1.1. The circular Jacobi ensemble. The theory of random unitary ma- 
trices was developed using the existence of a natural probability uniform 
measure on compact Lie groups, namely the Haar measure. The statisti- 
cal properties of the eigenvalues as well as the characteristic polynomial of 
these random matrices have played a crucial role both in physics (see [27J for 
an historical account) and in analytic number theory to model L-functions 
(see |20] and [21] where Keating and Snaith predict moments of L-functions 
on the critical line using knowledge on the moments of the characteristic 
polynomial of random unitary matrices). 

The circular unitary ensemble (CUE) is U{n), the unitary group over 
C", equipped with its Haar probability measure fJ-u{n)- Weyl's integration 
formula allows one to average any (bounded measurable) function on U{n) 
which is conjugation-invariant 

//dM^(n) = ;^/---/|A(e'^S...,e''")P/(diag(e'^S...,e'^"))^...|^, 

(l.f) 
where A(e'^^, . . . ,e^^") = ni<i<fc<n(^'^'° ~ ^'^^) denotes the Vandermonde 
determinant. 

The circular orthogonal ensemble (COE) is the subset of U{n) consisting 
of symmetric matrices, i.e. U{n)/0{n) = {VV^;V G U{n)} equipped with 
the measure obtained by pushing forward iJijj{n) by the mapping V i->- VV^ . 
The integration formula is similar to (jl.ip but with |A(e'^^, . . . , e'^")p re- 
placed by |A(e^^, . . . , e ")| and with the normalizing constant changed ac- 
cordingly. 

For the circular symplectic ensemble (CSE), which will not be recalled 
here, the integration formula uses |A(e^^^, . . . , e'^")|^. 

Dyson observed that the induced eigenvalue distributions correspond to 
the Gibbs distribution for the classical Coulomb gas on the circle at three 
different temperatures. More generally, n identically charged particles con- 
fined to move on the unit circle, each interacting with the others through 
the usual Coulomb potential — log \zi — Zj\, give rise to the Gibbs measure at 
temperature 1//3 (see the discussion and references in [2^ and in p^ chap. 
2): 

KU) = ctijf{e''\---.^'-)\^{e''\---.^'")fde^...den, (1.2) 

where CqL is a normalizing constant chosen so that 

4'2(ei,...,e„) = cg|A(e'^S...,e'^")|^ (1.3) 
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is a probability density on (0, In)"' and where / is any symmetric function. 
The unitary, orthogonal and syniplectic circular ensembles correspond to 
matrix models for the Coulomb gas at /3 = 1,2,4 respectively, but are there 
matrix models for general /? > for Dyson's circular eigenvalue statistics? 

Killip and Nenciu [24] provided matrix models for Dyson's circular en- 
semble, using the theory of orthogonal polynomials on the unit circle. In 
particular, they obtained a sparse matrix model which is five-diagonal, called 
CMV (after the names of the authors Cantero, Moral, Velasquez |9j). In this 
framework, there is not a natural underlying measure such as the Haar mea- 
sure; the matrix ensemble is characterized by the laws of its elements. 

There is an analogue of Dyson's circular ensembles on the real line: the 
probability density function of the eigenvalues {xi, . . . ,Xn) for such ensem- 
bles with temperature 1//3 is proportional to 

n 

\A{x,,...,Xn)flle-^V' (1.4) 

i=i 
For (3 = 1,2 or 4, this corresponds to the classical Gaussian ensembles. 
Dimitriu and Edelman [H] gave a simple tridiagonal matrix model for (II. 4p . 
Killip and Nenciu [21], gave an analogue matrix model for the Jacobi mea- 
sure on [—2, 2]", which is up to a normalizing constant, 

n 
|A(xi, ...,Xn)f Y[{2 - Xj)''{2 + xj)''dxi . . . dXn, (1-5) 

i=i 
where a,b > 0, relying on the theory of orthogonal polynomials on the unit 
circle and its links with orthogonal polynomials on the segment. When o 
and b are strictly positive integers, the Jacobi measure (jl.Sp can be inter- 
preted as the potential |A(xi, . . . , Xn+a+b)\^ on [—2, 2]"+"+^ conditioned to 
have a elements located at 2 and b elements located at —2. Consequently, 
the Jacobi measure on the unit circle should be a two parameters extension 
of (|1.3p . corresponding to conditioning to have specific given eigenvalues. 
Such an analogue was defined as the circular Jacobi ensemble in ^^ and 
|15| . If (5 G M, we recover the cJUE as in Witte and Forrester |37| . 

(n) 

Definition 1.1. Throughout this paper, we note h^ L the probability density 
function on (0, 2tt) given by: 



with (5 e C, d\c{5) > -i. 



/ig(0i, . . . , 0.) = cg|A(e^^S . . . , e^'-)f n(l - e"''0'(l - ^'^f (1-6) 



If (5 € ^N, this measure coincides with ()1.3p conditioned to have eigen- 



2' 

values at 1. For /? = 2, such measures were first considered by Hua [19j and 
Pickrell [30], [H]. This case was also widely studied in [29] and [6] for its 
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connections with the theory of representations and in |3i for its analogies 
with the Ewens measures on permutation groups. 

One of our goals in this paper is to provide matrix models for the circular 
Jacobi ensemble, i.e. a distribution on U{n) such that the arguments of the 
eigenvalues (e^^^, . . . ,e'^") are distributed as in (|1.6p . One can guess that 
additional problems may appear because the distribution of the eigenvalues 
is not rotation invariant anymore. Nevertheless, some statistical informa- 
tion for the circular Jacobi ensemble can be obtained from Dyson's circular 
ensemble by a sampling (or a change of probability measure) with the help 
of the determinant. More precisely, let us first define the notion of sampling. 

Definition 1.2. Let (X, J-", /x) be a probability space, and h : X >-^ M"*" a 
measurable and integrable function with E^(/i) > 0. Then a measure fi' is 
said to be the /i-sampling of /i if for all bounded measurable functions / 






Em'(/) 



If we consider a matrix model for HqI , we can definqj a matrix model for 

^5 ^y ^^^ means of a sampling, noticing that when the charges are actually 
the eigenvalues of a matrix U, then (|1.3|) differs from (|1.6|) by a factor which 
is a function of det(Id — U). We define det^ for a unitary matrix U as 

detsiU) = det(Id - C/)W(Id - 17)^ , 

and we will use this det^-sampling. 

Actually we look for an effective construction of a random matrix, for 
instance starting from a reduced number of independent random variables 
with known distributions. Notice that in the particular case /3 = 2, the 
density /io,2 corresponds to eigenvalues of a matrix under the Haar measure 
on U{n) and the det^-sampling of this measure is the Hua-Pickrell measure 
studied in our previous paper [8]. 

1.2. Orthogonal polynomials on the unit circle. We now wish to out- 
line the main ingredients which are needed from the theory of orthogonal 
polynomials on the unit circle to construct matrix models for the general 
Dyson's circular ensemble. The reader can refer to ^35j and |36J for more 
results and references; in particular, all the results about orthogonal poly- 
nomials on the unit circle (named hereafter OPUC) can be found in these 
volumes. 

Let us explain why OPUC play a prominent role in these constructions. 
In all this paper, D denotes the open unit disk {z € C : |2;| < 1} and T the 
unit circle {z € C : \z\ = 1}. Let {7i,u,e) be a triple where Ti is a Hilbert 
space, u a unitary operator and e a cyclic unit vector, i.e. {u^ e}'^_^ is 
total in T-L. We say that two triples (■H,n,e) and {IC,v,e') are equivalent if 
and only if there exits an isometry k : % ^s- K, such that v = kuk^^ and 



for 5He(<5) > —1/2, due to an integrability constraint 
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e' = ke. The spectral theorem says that for each equivalence class, there 
exists a unique probability measure /x on T such that 

{e,u''e)n= z^dfi{z) , k = 0,±l,.... 

Ji 

Conversely, such a probability measure n gives rise to a triple consisting of 
the Hilbert space L^(T, ;u), the operator of multiplication by z, i.e. h i— ?• 
{z I— 7- zh{z)) and the vector 1, i.e. the constant function 1. When the space 
% is fixed, the probability measure ^ associated with the triple {7i,u,e) is 
called the spectral measure of the pair (n, e). 

Let us consider the finite n-dimensional case. Assume that u is unitary 
and e is cyclic. It is classical that u has n different eigenvalues (e'^-'',j = 
l,...,n) with 6j G [0,2tt). In any orthonormal basis whose first vector 
is e, say (ei = e, 62, . . . ,e„), u is represented by a matrix U and there is 
a unitary matrix 11 diagonalizing U. It is then straightforward that the 
spectral measure is 



n 



where the weights are defined as ttj = \{ei,Iiej)\'^. Note that ttj > because 
a cyclic vector cannot be orthogonal to any eigenvector (and we also have 
J21=i ^j = 1 because 11 is unitary). The eigenvalues (e'^^', j = 1, . . . n) and 
the vector (tti, . . . ,7r„) can then be used as coordinates for the probability 
measure /x. 

Keeping in mind our purpose, we see that the construction of a matrix 
model from a vector (e^^^, j = 1, . . . ,n) may be achieved in two steps: first 
give a vector of weights (vri, . . . , -Kn), then find a matricial representative of 
the equivalence class with a rather simple form. The key tool for the second 
task is the sequence of orthogonal polynomials associated with the mea- 
sure /i. In -L^(T, /i) equipped with the natural basis {1, z, z"^, . . . , z"^"^}, the 
Gram-Schmidt procedure provides the family of monic orthogonal polyno- 
mials $0) • • • ) ^n-i- We can still define $„ as the unique monic polynomial 
of degree n with || $„ \\l'2{t,ij.)= 0, namely 

n 

(1.8) 



(1.9) 



(1.10) 

The coefficients a^-'s (0 < j < n — 1) are called Verblunsky coefficients and 
satisfy the condition ao, ■ ■ ■ , an-2 G ^ and Qn-i € T. 









The $fc's {k = 


= 0,.. 


. , n) obey the Szego recursion relation: 

<^j+i{z) = z<i>j{z) - aj<^*{z) 


where 








<!>*{z) = z^<!>,{-z'^). 
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When the measure /x has infinite support, one can define the family of 
orthogonal polynomials {^n)n>o associated with /x for all n. Then there are 
infinitely many Verblunsky coefficients {on) which all lie in B. 

Verblunsky's Theorem (see for example |35j . |36j ) states that there is a 
bijection between probability measures on the unit circle and sequences of 
Verblunsky coefficients. The matrix of the multiplication by z in L^(T, ^), 
in the basis of orthonormal polynomials, has received much attention. This 
unitary matrix, noted G{ao, • • • , On-i), called GGT by B. Simon (see Chap- 
ter 4.1 in [35] for more details and for an historical account), is in the 
Hessenberg form: all entries below the subdiagonal are zero, whereas the 
entries above the subdiagonal are nonzero and the subdiagonal is nonnega- 
tive (see formulae (4.1.5) and (4.1.6) in [35] for an explicit expression for the 
entries in terms of the Verblunsky coefficients, or formula (jl.lip in Lemma 
O below). 

For H a. n X n complex matrix, the subscript Hij stands for {ei,H{ej)), 
where {x,y) = J^k^i'^kUk- Killip and Nenciu state that any unitary ma- 
trix in the Hessenberg form with nonnegative subdiagonal is the matrix of 
multiplication by z in L^(T, /i) for some measure /i. More precisely, from 
Killip-Nenciu [23], Lemma 3.2, we have 



Lemma 1.3. Let fi be a probability measure on T supported at n points. 
Then the matrix H of f{z) \-^ zf{z) in the basis of orthonormal polynomials 
of Lp'iT,^), is in the Hessenberg form. More precisely 

-tti^iOj U^pl] pp ifi<j + l 

Pj-i ifi=j + l, (1.11) 

ifi>j + l 



Hi+ij+i 



with pj = Y^l — |ajp and a_i = — 1, the au's being the Verblunsky coeffi- 
cients associated to p.. Conversely, if uq, . . . ,0^-1 are given in D^~^ x T, 
and if we define the matrix H by \1.11\) . then the spectral measure of the 
pair {H,ei) is the measure p whose Verblunsky coefficients are precisely 
OiQ, . . . ,an-i- 

Besides, there is a very useful decomposition of these matrices into prod- 
uct of block matrices, called the AGR decomposition by Simon ([33])) after 
the paper [3]. For < k < n — 2, let 



e(^-)(a) = idfcef "'^ P'' ) (Bid 

\ Pk -otk J 



n-k-2- 



and set 0^" ^'{an-i) = Idn-i ® (an-i), with |a„_i[ = 1. Then the AGR 
decomposition states that ([33] Theorem 10.1) 

g{ao, . . . ,a„-i) = e(o)(ao)e(i)(ai) . . . e("-i)(a„_i) . 

Now we state a crucial result of Killip and Nenciu which enabled them to 
obtain a matrix model in the Hessenberg form for Dyson's circular ensemble. 
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The challenge consists in randomizing (oq, • • • , On-i) in B"^"'^ x T in such a 
way that the angles {6i, . . . , On) of the spectral measure 

n 

i=i 

have the density /iq i (see (|1.6p ). To make the statement precise, let us 
introduce three definitions (for the properties of the Beta and Dirichlet dis- 
tributions, see 



Definition 1.4. For ai, 02 > 0, let Beta(ai, 02) be the distribution on [0, 1] 
with density 

r(ai+a2) ai^l.-l _^xa2-l 

Its generalization is the following. 

Definition 1.5. For n > 2 and oi, . . . , a„ > 0, let Dir(ai, . . . , a„) be the 
distribution on the simplex {(xi,...,x„) € [0,1]" : Y17=i^i — 1} with 
density 



r(ai + ... + a„) -pj- ^^k-i 
k=l 



r(ai)...r(a,^ Ai '^ 

If ai = . . . = a„, = a, it is called the Dirichlet distribution of order n > 2 
with parameter a > and denoted by Dirn(a). (For a = 1, this is the 
uniform distribution). 

Definition 1.6. For s > 1 let Vs be the probability measure on D with 
density 

i^(l_[^|2)(.-3)/2_ 

2tt 
It is the law of Re'^ where R and ip are independent, Tp is uniformly dis- 
tributed on (0, 2ir) and R^ has the Beta(l, (s — l)/2) distribution. We adopt 
the convention that vi is the uniform distribution on the unit circle. We 
denote by rj^l the distribution on B"^"*^ x T given by 

??i"j = 0^-oV^(„_fc_i)+i . (1.12) 

Proposition 1.7 (Killip-Nenciu [24J, Proposition 4.2). The following for- 
mulae express the same measure on the manifold of probability distributions 
on T supported at n points: 

r)l — n "■ 

-—\A{e''\. . . , e''-)f rr n^/^-'dOi . . . dOndni . . . d^n-i 

in the (9, vr) coordinates and 

n-2 , , 

JJ(1 _ |«^|2)(/3/2)(n-fc-l)-1^2^^ _ _ _ ^'«n-2|^ (1-13) 

fc=0 
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in terms of the Verblunsky coefficients. 

Proposition [TT71 may be restated as follows: to pick at random a measure // 
such that {ao, ■ ■ ■ , On-i) is rj^l distributed is equivalent to pick the support 

{9i, . . . ,9n) according to h^l (see (|1.3p ) and to pick the weights (vri, . . . , vr^) 
independently according to Dir„(/3/2). 

As a consequence, if one takes independent coefficients {ao, . . . , a„_i) such 
that Ok is i^i3(n-k-i)+i distributed for < A; < n — 1, then the GGT matrix 
G{aQ, . . . , On-i) will be a matrix model for Dyson's circular ensemble with 
temperature 1/(3 (see also Proposition 2.11 in [13j). Actually in |24j . Killip 
and Nenciu provide a matrix model which is much sparser (five-diagonal) as 
shall be explained in Section [H 

Let us now define the laws on U{n) which we will consider in the sequel. 

(n) 

Definition 1.8. We denote by CJq i the probability distribution supported 
by the set of n x n GGT matrices of the form (jl.lip . corresponding to the 
law of Q{ao, . . . ,0^-1) defined above. We denote by Cjf^l the probability 
distribution on U{n) which is the det^-sampling of CJq i. 

The standard GGT approach is not sufficient to produce matrix mod- 
els for the circular Jacobi ensemble because, as we shall see in Section [3l 

(n) 

under the measure CJ^J, the Verblunsky coefficients are not independent 
anymore. To overcome this difficulty, we associate to a measure on the unit 
circle, or equivalently to its Verblunsky coefficients, a new sequence of co- 
efficients (7/c)o</c<n-i) which we shall call deformed Verblunsky coefficients. 
There is a simple bijection between the original sequence (afc)o<fc<n-i and 
the new one (7fc)o<fc<n-i- These coefficients satisfy among other nice prop- 
erties that \ak\ = |7fc|, and that they are independent under CJ^o (and for 
6 = the a^s and the 7fc's have the same distribution). They have a geo- 
metric interpretation in terms of reflections: this leads to a decomposition 
of the GGT matrix Q{ao, . . . On-i) as a product of independent elementary 
reflections (constructed from the 7fc's). The explicit expression of the den- 
sities allows an asymptotic study (as n —^ 00) of the 7a;'s, and consequently 
of the spectral measure, and flnally of the empirical spectral distribution. 

1.3. Organization of the paper. In Section 2, after recalling basic facts 
about the reflections introduced in [8], we deflne the deformed Verblunsky 
coefficients (7fc)o<fc<n-i and give some of its basic properties. In particular 
we prove that the GGT matrix G{ao, . . .On-i) can be decomposed into a 
product of elementary complex refiections (Theorem 12. 8p . 

In Section 3, we derive the law of the 7fc's under CJ^"^ (Thorem l3.2p : in 
particular we show that they are independent and that the actual Verblunsky 
coefficients are dependent if (5 7^ 0. We then prove an analogue of the above 
Proposition 11.71 on the {6,tt) coordinates of /Li (Theorem 13. 3p . 
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In Section 4, we propose our matrix model (Theorem 14. ip . It is a modifi- 
cation of the AGR factorization, where we transform the Ofc's so that they 
become reflections : 

E(''\a) = idfc e ( ^ _f^^ ) e id„_,_2, 

with e^''' = jE^- Of course the CMV representation [9], which is five- 
diagonal, is also available, but this time the a^'s are not independent. Using 
the following elementary fact proven in Section 2, 

n-l 

$„(l) = det(Id-?7) = [](l-7fc), 

fc=0 

we are able to generalize our previous results in [7] and [8] about the de- 
composition of the characteristic polynomial evaluated at 1 as a product of 
independent complex variables (Proposition I4.3( ) . 

In Section 5, we study asymptotic properties of our model as n — > cxd, 
when 5 = /3nd/2, with IHed > 0. We first prove that the Verblunsky coeffi- 
cients have deterministic limits in probability. This entails that the spectral 
measure converges weakly in probability to the same deterministic measure 
(denoted by /x^) which is supported by an arc of the unit circle (Theorem 
15. ip . Besides, we consider the empirical spectral distribution (ESD), where 
the Dirac masses have the same weight 1/n. Bounding the distances be- 
tween both random measures, we prove that the ESD has the same limit 
(Theorem 15. 4p . Moreover, starting from the explicit joint distribution ()1.6p . 
we prove also that the ESD satisfies the large deviation principle at scale 
{(i/2)n? whose rate function reaches its minimum at /Xd (Theorem 15. 5p . 

2. Deformed Verblunsky coefficients and reflections 

In this section, we introduce the deformed Verblunsky coefficients and we 
establish some of their relevant properties, in particular a geometric interpre- 
tation in terms of reflections. One remarkable property of the Verblunsky co- 
efficients, as it appears in Proposition !!.?! is that they are independent under 
CJq o. As we shall see in Section[3l this does not hold anymore under CJ^'o. 
This motivated us to introduce a new set of coefficients, (70, . . . , 7„_2, 7n-i), 
called deformed Verblunsky coefficients, which are uniquely associated with 
a set of Verblunsky coefficients. In particular, 7^ € D for < fc < n — 2, 
7n_i G T and the map (70, . . . , 7n-i) ^^ (aO) • • • i ctn-i) is a bijection. More- 
over, the characteristic polynomial at 1 can be expressed simply in terms of 
(7o,...,7„_i). 

2.1. Analytical properties. Let ^ be a probability measure on the unit 
circle supported at n points. Keeping the notations of the introduction, we 
let (^A:(-2;))o<fc<n denote the monic orthogonal polynomials associated with 
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H and {ak)o<k<n-i its corresponding set of Verblunsky coefficients through 
Szego's recursion formula (II. 9p . The functions 

bk{z) = ^y k<n-l (2.1) 

are known as the inverse Schur iterates ([36j p. 476, after Khrushchev [23] 
p. 273). They are analytic in a neighborhood of D and meromorphic in C. 
Each bk is a finite Blashke product 

n 

h{z) = n 



1 



where zi, . . . ,Zk are the zeros of $fc. Let us now explain the term "inverse 
Schur iterate". 

The Schur function is a fundamental object in the study of the orthogonal 
polynomials on the unit circle. Let us briefly recall its definition (see [35] 
or [M] for more details and proofs): if // is a probability measure on the 
unit circle (supported at finitely many points or not), its Schur function 
/ : D — > D is defined as: 

f{z) = -^^^ where F{z) = /^d/x(e^^). (2.2) 



z F{z) + 1 ^ ^ J e''^ - z 

It is a bijection between the set of probability measures on the unit circle and 
analytic functions mapping B to B. The Schur algorithm (which is described 
in ^35| or [34] p. 438) allows to parametrize the Schur function / by a sequence 
of so-called Schur parameters, which are actually the Verblunsky coefficients 
associated to /i (Geronimus theorem). In particular, there are finitely many 
Verblunsky coefficients (or equivalently the measure /i is supported at n 
points) if and only if / is a finite Blaschke product. The name "inverse 
Schur iterate" ([22j) for bk comes from the result (1.12) of the latter paper 
where bk is identified as the Schur function corresponding to the "reversed 
sequence" {—ak-i-, • • • , — ao> 1) (see also [36] Prop. 9.2.3). 

Let us define our sequence of functions, which shall lead us the deformed 
coefficients. 

Definition 2.1. If ^ is supported at n points and with the notation above, 
define ^k[z) for < A; < n — 1, as: 

<^k+i{z) .„„. 

7.(^) = ^-^^(^. (2.3) 

From the Szego's recursion formula (II. DP and notation (12. ip . this is equiva- 
lent to 

lk{z) = -^ , (2.4) 

bk{z) 

so that 7yfc is meromorphic, with poles in B and zeros lying outside B. 
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The next proposition shows how the functions 7fc(-?) can be defined recur- 
sively with the help of the coefficients at- As a consequence, we shall see 
that the 7^(2) are very closely related to a fundamental object in the theory 
of random matrices: the characteristic polynomial. 

Proposition 2.2. For any 2; € C, 70 (-z) = ao O'l^-d the following decomposi- 
tion for ^ki^) holds: 

fc-i 
^k{z) = ll{z-j,{z)) , k = l,...,n. (2.5) 

j=0 
The 7fc(z) 's may be also defined by means of the a's through the recursion : 



7*(3) = 7t(J-')^ (2^7) 

Proof. The first claim is an immediate consequence of (j2.3p . Now, using 
^k{z) = W]zl{z - lj{z)), we obtain 

fe-i 
K{z) = \{{l-z^,{z)), 

j=0 

and hence (we use (|2.4p ) 

1- z^jjz) 



lk{z) = akT\ . 



D 



Note that when \z\ = 1, J7a;(-z)| = |afc|. Combined with the above propo- 
sition, this leads us to introduce the following set of coefficients. 

Definition 2.3. Define the coefficients (7fc)o<fc<n-i by 

7fc:=7fc(l), k = 0,...,n-l. (2.8) 

We shall refer to the 7fc's as the deformed Verblunsky coefficients. 

Proposition 2.4. The following properties hold for the deformed Verblunsky 
coefficients: 

a) For all < k < n — 1, {j^l = \oik\, <ind in particular jn-i G T; 

b) 7o = ao and 

fe-i 

7fc = afce^^'=-i , e'^''-- = T\—^ , (fe = 1, . . . ,n - 1) . (2.9) 

i — 7o 

The last term is special. Since \an^i\ = 1, we set Un-i = e^"^"'^ , so that 
7„_i = e'(-'^"-i+^"-2) := e'^"-i . (2.10) 
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c) Let fi be the spectral measure associated to {U,ei), U G U{n). Then 
^n{z) is the characteristic polynomial ofU, in particular, 

n-\ 

$„(l) = det(Id-[/) = J];(l-7fc). (2.11) 

fc=0 

Proof. All the results are direct consequences of the definition 12.31 and the 
formulae in Proposition 12.21 evaluated at 1. D 

Remark. In |26j . Killip and Stoiciu have already considered variables which 
are the complex conjugate of our deformed Verblunsky coefficients as aux- 
iliary variables in the study of the Priifer phase (Lemma 2.1 in [26]). Nev- 
ertheless, the way we define them as well as the use we make of them are 
different. 



Remark. The formula (12. 9p shows that the jkS can be obtained from the 
ak's recursively. Hence starting from a spectral measure associated to a 
unitary matrix, one can associate with it the Verblunsky coefficients and 
then the deformed Verblunsky coefficients. Conversely, one can translate 
any property of the deformed ones into properties for the spectral measure 
associated with it by inverting the transformations ()2.9p . 

Remark. The distribution of the characteristic polynomial of random unitary 
matrices evaluated at 1, through its Mellin- Fourier transform, plays a key 
role in the theory of random matrices, especially through its links with 
analytic number theory (see [28] for an account). In [7] it is proven that it 
can be decomposed in law into a product of independent random variables 
when working on the unitary and orthogonal groups endowed with the Haar 
measure; since we will prove in Section 3 that the 7fc's are independent under 

(n) 

CJgl, then we can conclude that this latter result holds for any circular 
Jacobi ensemble. 

2.2. Geometric interpretation. We give a connection between the coef- 
ficients (7fc)o<fc<n-i and reflections defined just below. This allows us to 
obtain a new decomposition of the GGT matrix associated with a measure 
fj, supported at n points on the unit circle as a product of n elementary 
reflections. 

Many distinct definitions of reflections on the unitary group exist, the 
most well-known may be the Householder reflections. The transformations 
which will be relevant to us are the following ones. 

Definition 2.5. An element r in U{n) will be referred to as a reflection if 
r — Id has rank or 1. 

If w G C", we denote by {v\ the linear form w i— )■ {v,w). The reflections 
can also be described in the following way. If e and m ^ e are unit vectors 
of C", there is a unique reflection r such that r(e) = m, and 

r = Id-- r(m-e)((m-e)| . (2.12) 

-L \ III/ 1 t- / 
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Let F := span{e, m} be the 2-diniensional vector space which is spanned by 
the vectors e and m. It is clear that the reflection given by formula ()2.12p 
leaves F invariant. Now set 

7 = (e,m) , p = yr^bF , e'^ = ^ , (2-13) 

1-7 

and let g G F be the unit vector orthogonal to e obtained by the Gram- 
Schmidt procedure. Then in the basis {e,g) of F, the matrix of the restric- 
tion of r is 

=(7) - ( I i^ ) . (2.14) 

Conversely, for 7 € B, such a matrix represents the unique reflection in C^ 
provided with its canonical basis, mapping ei onto 761 -|- y^l — |7pe2. The 
eigenvalues of r are 1 and — e"^. 

Let n be a unitary operator in C" and e a unit cyclic vector for u. We 
define n reflections ri, . . . ,r„ recursively as follows. Let (ei, . . . ,en) be the 
orthonormal basis obtained from the Gram-Schmidt procedure applied to 
(e, ue, . . . ,it"'~^e). 

Let ri be the reflection, mapping e = £1 onto ue = u£i. More generally, 
for /c > 2 let r^ be the reflection mapping e^ onto f^^-^rj^^^ ■ ■ -^i ue^. We 
will identify these reflections and establish the decomposition of u. Following 
the basics recalled about GGT matrices in the introduction, we note that the 
matrix of u in the basis (ei , . . . , e„) is the GGT matrix associated to the mea- 
sure fi, i.e. the matrix Q{ao,--- ,a„-2,an-i)) where {ao,-" ,Oin-2,Oin~i) 
are the Verblunsky coefficients associated with the measure /x. We will use 
formula (4.1.6) of [35] or formula (jl.lip of (our) Lemma 11.31 for the identi- 
fication of scalar products. 

Proposition 2.6. (1) For every 1 < k < n — 1, the reflection r^ leaves 

invariant the n— 2 -dimensional space Spanjei, . . . ,ek-i,£k+2, ■ ■ ■ ,£n}- 
The reflection rn leaves invariant Span{£i, . . . ,e„_i}. 
(2) The following decomposition holds : 

u = ri---rn. (2-15) 

Proof. (1) In view of Section 2.1, it is enough to prove that for j ^ {k, k + 1}, 
the vectors £j and rkSk are orthogonal. 

For k = 1, {£j,ri£i) = {ej, usi) = as soon as j > 3 from (II. IIP . 

Assume that for every ^ < /c — 1, the reflection ri leaves invariant 
Spanjei, . . . , q_i, Q+2, • • • , £n}- For every j = 1, . . . , n, we have 

{ej,rk£k) = {£j,r'i:\rj:^^...r^'^u£k) = {ri ■ ■ -rk-iSj^uek) ■ (2.16) 

For j > k + 2, hy assumption, the reflections ri, . . . , r^-i leave invariant 
Ej, SO that the above expression reduces to {Ej,uEk) which is again by 

dni]). 

For j = k — 1, we have n • • • rk-iEk-i = uEk~i by deflnition of r^-i, so 
that (I2.16P gives {Ek-i,rkEk) = {uEk-i,uEk), which is since u is unitary. 
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For j < k — 1, by assumption, the reflections Vj^i, . . . , r^-i leave invariant 
Ej, SO that the right hand side of (|2.16p reduces to {ri ■ ■ -rjej^uEk)- By 
definition of rj, it is {uej,uek) which is 0. 

(2) For k fixed, it is clear from (1) that ri • • • r„ejt = ri ■ ■ ■ rk£k which is 
uEk by definition of r^. D 

Proposition 2.7. For k = 1, . . . ,n — 1, the matrix of the restriction of rj. 
to the basis {ek,£:k+i) is H(7fc_i) as defined in {2.14\) - In particular 

{£k,rk£k) =7fc-i- (2-17) 

The restriction of rn to Ce„ is the multiplication by 7n-i- 

Proof. Note that for every k < n — 1 

(efc+i, r^efc) = {n ■ ■ ■ rk-iek+i,uek) = {ek+i,uek) = Pk~i ■ (2.18) 

Since r^ is a reflection acting on the subspace Span{ek,£k+i}i identities 
(j2.18p and (j2.17p entail that the matrix representing r^ in the basis (ei, . . . , e^) 
is precisely E{'~fk-i) (see (|2.14p ). It is then enough to prove ()2.17p . 
For /c = 1 it is immediate that: 

(ei,riei) = (ei,uei) = ap = 7o • 

Let us proceed by induction. For j > 1 set qj := {ej,rj£j). Assume that qj = 

7i-i for j < k We have g^+i = {£k+i,rk+i£k+i) = (n . . . ^fc^fc+i, ^i^fc+i). 
Equation (j2.12p implies 

I 

rk£k+i = efc+i - :; = {rkek - ek){rk£k,ek+i), (2.19) 

1 - Ik-i 

and since rj££ = e^ for £ > j + 2, we get: 

ri.-.rkEk+i = Ek+i - z (n ---rkek - ri . . .rk-iek){uek,£k+i)- 

J- - 7fc-i 



Now, it is known that {uek,£k+i) = {£k+i,uek) = Pk- If we set t^i = £i, 

Vj =ri... rj^iEj , aj = ^~_^ , Wj+i = Ej+i - GjUEj 

1 - 7i-i 

we get the recursion 

Vj+i = UjVj + Wj+i , {j < k) , (2.20) 

which we solve in : 

k fc+l k 

i=i e=2 j=e 

Taking the scalar product with uEk+i yields 

k fc+l k 

Qk+1 = ( JJaij(ei,tiefc+i) + ^ [^Y{ajj{wi,uEk+i). 

j=l £=2 j=£ 
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But {we,uek+i) = {ee,uek+i) - ae-i{uee^i,uek+i), and since i < k + 1, we 

have 

fc-i 

{wi,uek+i) = {ei,uek+i) = -akae_2 JJ Pm, 

m=e-i 

which yields (with a_i = —1) 

Qk+l 



Oik 



fc+l k fc-1 

e=i j=e m=i-\ 

fc+l /c— 1 £—3 _ 

£=1 m=f-l j=0 Hs=OV^ is) 

fe+l fc-1 £-3 fc-1 £-2 

lls=0'>^ Is) i^i rn=i-2 j=Q m=i-l j=0 

W (1 - 7.) 

and eventually g^+i = 7fc- D 

Now, we can summarize the above results in the following theorem. 

Theorem 2.8. Letu S U{n) ande a cyclic vector for u. Let fi be the spectral 
measure of the pair {u,e), and (ao, • • • ,an-i) its Verblunsky coefficients. 
Let {ei, . . . ,£n) be the orthonormal basis obtained from the Gram-Schmidt 
procedure applied to {e,ue, . . . ,u"'~^e). Then, u can be decomposed as a 
product of n reflections {rk)i<k<n-' 

u = ri...rn (2.22) 

where ri is the reflection mapping £i onto uei and by induction for each 
2 < k < n, rk maps Ek onto r^_^r^_2 . . .r^ uEk- 

This decomposition can also be restated in terms of the GGT matrix : 

g{ao, ■■■ , an-i) = hW(7o)H(i)(7i) . . . E^^^^'H^n-i) , (2.23) 

where for < k < n — 2, the matrix H*- •* is given by 

H('=-i)(7fc_i) = Idfc_i © S(7fc„i) © Id„_fc_i, (2.24) 

with 2(7) defined in ^2.14^ . For k = n — 1, 

S("-I)(7„_i) = Id„_i © (7„_i). (2.25) 

3. Deformed Verblunsky coefficients and independence 

We now use the point of view of sampling (or change of probability mea- 
sure) to compute the distribution of the deformed Verblunsky coefficients 
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(n) 

under CJ^ i. Let us first remark that, if ttie a^'s are independent and witli 
rotational invariant distribution, then from (|2.9p 



(ao,...,an-i) = (7o,---,7n-i)- (3.1) 

This is the case under CJq^. 

We first prove that when 6^0 the Verblunsky coefficients are not in- 
dependent anymore by studying the simple case n = 2, /3 = 2, and then 
we compute the distribution of (70, . . . ,7n-i) under CJ^"i. We then show 
that under this distribution, the weights of the measure associated to the 
Verblunsky coefficients {oq, . . . ,an-i) are independent from the points at 
which the measure is supported and follow a Dirichlet distribution. 

Let (5 G C such that 9\e6 > -1/2. The formula 

^(,) ^ r(i + ^)r(i + ^) ^^ , 

defines a probability density with respect to the Haar measure on T, which 
is discontinuous at 1 when 3m 5 7^ (see [8]). 

When (3 = 2 and (5 7^ 0, the Verblunsky coefficients are dependent. Indeed, 
let M E U{2) with Verblunsky coefficients oq and ai. Then 

det(Id - M) = [1 - ao - ai(l - oq)], 

(2) 
with |ao| < 1 and [ai| = 1. Under CJq2, the variables ao and ai are 

independent and uniformly distributed on D and T respectively (see [23] or 

Proposition 1 1.7p . The CJ^ 2 is a det^ sampling of CJq 2 (see the Introduction 

for the definition and notation for the det^- sampling). So, the joint density 

of (a, if) on D X T is proportional to 

f{ao,ai) = [1 -ao -ai(l -ao)]''[l -ao -ai(l -ao)]"^ 
= (1 - a-o)^"(l - ao)'^[l - jaif[l - 701]^ . 

where 

1 — a 

7 = :; -■ 

1 — a 

It is then clear that the conditional density of ai given ao is 

ai h^ A('')(7a-i) , ai G T, 

and this last quantity does not depend on ao (i.e. on 7) if and only [{5 = 0. 
Otherwise the Verblunsky coefficients ao and ai are dependent. 

The next theorem illustrates our interest in the deformed Verblunsky co- 

(n) 

efficients: under CJ^ i, they are independent. For the proof of this theorem, 
we shall need the following lemma which will also be useful when we study 
limit theorems: 
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Lemma 3.1. Let s,t,£eC such that: D\t{s + ^ + 1) > 0, D\t{t + £ + !)> 0. 
Then, the following identity holds: 

Proof. A Taylor expansion yields 

^-^ n\vn\ 

m,n>0 

with z = pe'^ and < p < 1. We obtain by integrating 

(1 - \zfY-\l - zYil - zfd^Z = 2vry izfMz^k [\l _ ^2)£-1^2n+l^^ 
n>0 ^ ' 

where 2^1 is the classical hyper geometric function (see [4]) and an applica- 
tion of Gauss formula (see [4j) shows that the last expression is exactly the 
right hand side of (j3.3p . D 

Theorem 3.2. Let 5 e C with «He 5 > -1/2 and /3 > 0. Set /?' = /3/2. 
Under CJf'l, the distribution of (701 • • • i7n.-i); denoted hereafter rj^I, is 
the following: 

(1) the variables 70, • • • ,7n-2i7n-i o,re independent ; 

(2) for k = 0, . . . , n — 2 the density of 'jk with respect to the Lehesgue 
measure d^z on C is 

Cfc,„(5) (1 - |^|2)/3'("-'=-l)-l (1 _ ^)5(i _ ^)<5i^(^) ^ 

where 

_ T{P'{n - fc - 1) + 1 + 5)T{p'{n - k - I) + I +-5) _ 
""^'""^ ' ^r(/3'(n -k- l))r(/3'(n -k -I) + 1 + 5 + 5) ' ^"^ 

('^j the density of 7„_i mi/i respect to the Haar measure on T is \^ ' 
(given by [37^) ). 

(n) 

Proof. The distribution of the a's in the /3-circular unitary ensemble is Vq s- 
More precisely, as seen in Definition 11.51 thev are independent and if 

ak = Rke'^" {0<k<n-2) , q„_i = e'^-^ , 

then Rfc and ipk are independent, ipk is uniformly distributed and r| has the 
Beta(l,/3'(n — A; — 1)) distribution. Moreover On-i is uniformly distributed 
on T. From (|2.11|) . the sampling factor is 

n-2 

det(Id - C/)^"det(Id - UY = (1 - 7n-i)'(l - In-if H^^ " ^'^)'(1 " ^'^)' • 

k=0 



18 P. BOURGADE, A. NIKEGHBALL AND A. ROUAULT 

SO that, under CJj^L the density of (Rq, . . . , Rn-2-,ipo^ ■ ■ ■ ■, V'n-i) is propor- 
tional to 

n-2 

A(^) (7n-i) n (1 - 4/^""-'-'^-'^^ - 7fc)'(l - %)'tio,i){Rk) , 

k=0 

with 

7fc = Rfcc'^^ (0 < fc < n - 2) , 7„_i = e'^"-i . 
Thanks to the relations (j2.9p and (|2.10p . the Jacobian matrix of the mapping 

(RO, . . . , Rn-2, -00, • • • , ■0n-l) ^ (R-O, • • • , Rn-2, ^0, • • • , ^n-l) 

(n) 

is lower triangular with diagonal elements ±1, so that, under CJ^i, the 
density of (Rq, . . . , Rn-2i ^O) • • • i ^n-i), is proportional to 

n-2 

A(')(7n-i) 11(1 - Ri)/^'("-i-'=)-iR,(l - 7fc)'"(l - 7fc)'l(o,i)(Rfc) , (3.5) 

fc=0 

which proves the required independence and the expression of the distribu- 
tions, up to the determination of Cfc^„((^). This quantity is obtained by taking 
i = j3'{n — k — 1) ^ s = 5 ,t = 5 in (j3.3p . which gives (j3.4p and completes the 
proof of the Theorem. D 

Starting with a set of deformed Verblunsky coefficients, with distribution 
rjf'a, we obtain the coefficients (qq, • • • , an-i) by inverting formula ()2.1ip . 
These are the coordinates of some probability measure // supported at n 
points on the unit circle: 

n 
fc=l 

with TTfc > and Yl'k=i''^k = 1- The next theorem gives the distribution 
induced on the vector (vri, . . . , vr^, 6*1, . . . , 6l„) by (70, . . . , 7n-i)- 

Theorem 3.3. The following formulae express the same measure on the 
manifold of probability distribution on T supported at n points: 



K^^\A{e''\. . .,e''-)f ll{l-e~'''')'{l-e''^yi[4'-'d9, . . . d^.dvn . . . d7r._i 



fe=i fe=i 

in the {9,tt) coordinates and 

n—2 n— 1 

K^^ 11(1 - h,\Y'^-'-''>-' 11(1 - 7.)^"(1 - 7.)'Ao . . . d^n-2# 
fc=0 fc=0 

in terms of the deformed Verblunsky coefficients, with 7^-1 = e^"^. Here, 
K^g is a constant: 



<5,/3 2"-i7rr 1 + (5 + 5 ,■'■-'- '^ ' 
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with Cfc^„((5) given in Theorem 13.^1 Consequently, if {jq, . . . ,jn-i) is ry^l 
distributed, then (vri, . . . , 7r„) and {9i, . . . ,6n) are independent; the vector of 
weights (tti, . . . , 7r„) follows the Dirn{f3') distribution and the vector {9i, . . . ,9n 

has the density h^ I . 

Proof. In the course of this proof, we shah adopt the following point of view. 
Starting with a measure supported at n points on the unit circle, we associate 
with it its Verblunsky coefficients (ao) ■ ■ ■ > On-i) and then the corresponding 
GGT matrix which we note G for simplicity. Then ei is a cyclic vector for 
G and // is the spectral measure of (G, ei). Conversely, starting with the 
set of deformed Verblunsky coefficients with r]s^i3 distribution, we construct 
the coefficients {oq, . . . , On-i) with the transformations ()2.9p . then the GGT 
matrix associated with it and finally /x the spectral measure associated with 
this matrix and ei. 

We use the following well-known identity (see [35] or [24] Lemma 4.1): 

n n—2 

|A(ei^S...,ei^")|2J]^fe = n(l-|«fcl'r"'"'- (3-6) 

fc=l A:=0 

Since \^k\ = lo^fcl) we can also write 

n n—2 

|A(ei^^...,ei^")|2J]^, = ^(l-|7fcP^"'"'• (3-7) 

fc=l fc=0 

Moreover, from (12.111) . 

n n—1 

det(Id-G) = Hil-e''^) = n(l-7fc)- (3-8) 

k=l k=0 

In our setting, ttj is modulus squared of the first component of the i-th 
eigenvector of the matrix G. Now, define 

Qk = ^k, ioT k = l,...,n. 
It is known (see for example Forrester |13j , Chapter 2 and |14] Theorem 2) 
that the Jacobian of the map (oq, • • • , an-i) ^ (6*1, • • • On, Qi, ■ ■ ■ , Qn-i) is 
given by 

mZ'oil - \ak\') 

nn 
k=l1k 

Moreover, the map (70, . . . ,7n-i) '-> {oiq,- ■ ■ ,oin-i) is invertible and its 

Jacobian is 1, as already seen. The result now follows from simple integral 

manipulations combined with the identities (13. 7p and ()3.8p . D 



4. Matrix models for the circular Jacobi ensemble 

The results of the previous sections can now be used to propose some 
simple matrix models for the circular Jacobi ensemble. There are mainly 
two ways to generate matrix models for a given spectral measure encoded 
by its Verblunsky coefficients. 
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The AGR decomposition : iiU = eW(ao)6^^Hai) • • • e("-i)(a„_i) (the 
Ofc's are defined in the introduction), the Verblunsky coefficients for the 
spectral measure associated to (C/, ei) are precisely (ao) • • • i ctn-i) (see [3] 
or [33j Section 10). Therefore, taking independent a^'s with law r/^ i, the 
density of the eigenvalues of 

U = e(°)(ao)eW(ai) . . . G("-i)(a„_i) 

is proportional to |A(e^^, . . . ,e ")|'^. The matrix U obtained above is the 
GGT matrix associated with the a^'s. It is in the Hessenberg form. 
The CM V form : set 

C = e(o)(ao)0^^Ha2)... 
M = eW(ai)e(3)(a3)... 

Cantero, Moral, and Velazquez [9] proved that the Verblunsky coefficients 
associated to (>C7W,ei) are precisely (ao, • ■ • , an-i)- Therefore, taking as 

In) 

previously independent a^'s with distribution "% L the density of the eigen- 
values of the spectral law of £A^ is proportional to |A(e^^i, . . . , e'^")|^ (}24j). 
This matrix model is very sparse: it is pentadiagonal. 

We now propose a matrix model for the circular Jacobi ensemble: it is 
reminiscent of the AGR factorization with the noticeable difference that it 
is based on the deformed Verblunsky coefficients and actual reflections as 
defined in Section 2. 

Theorem 4.1. // (70, . . . ,7n-i) is Vsr distributed^ then with the notation 
of (2^ and IKMl . 

s(o)(7o)H«(7i)...h("-i)(7„_i) 

is a matrix model for the circular Jacobi ensemble, i.e. the density of the 
eigenvalues is h^l (see lil.6\) ). 

Proof. We know from Theorem 12.81 that 

g{ao, • • • , a„_i) = h(o)(7o)h(i)(7i) . . . h("-i)(7„„i). (4.1) 

We also proved in Theorem 13.31 that the set of deformed Verblunsky coeffi- 

(n) 

cients with probability distribution rj^ I induces a distribution on the eigen- 
values of the GGT matrix Q{aQ, ■ ■ ■ ,q„_i) which has exactly the density 
hs^l3. This completes the proof of the Theorem. D 

Remark. We now say a few extra words on the CMV form obtained by Killip 
and Nenciu in [24]. Cantero, Moral and Velazquez [9] introduced the basis 
XO) • • • ) Xn-i obtained by orthogonalizing the sequence 1, z, z~^, .... They 
prove that in this basis the matrix is pentadiagonal. We name this matrix 
C{ao, • • • , a„_i). It turns out that there exists a unitary P such that: 

Pg{ao, ■■■ , an-i)P* = C(qo, • • • , an-i) , Pfo = Xo- 
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The two pairs {G{ao,--- ,an~i),^o) and {C{ao,--- ,a„_i),xo) are equiv- 
alent, they admit the a^s as Verblunsky coefficients, and have the same 
spectral measure. We conclude that if we start with the 7fc's distributed 
as r)gl, and build the a^s by inverting the transformation (|2.9p . then 
C(ao,--- ,an-i) will be a matrix model for the circular Jacobi ensemble. 
But we do not know how to construct the CMV matrix from the 7fc's directly. 
We saw at the beginning of this section that Cantero et al. introduced the 
matrices C and Ai, as direct product of small blocks @^^\ak) and obtained 
C as C = C^A. It would be interesting to have an analogue construction 
based on the independent 7fc's. 

Theorem 12.81 which is a deterministic result, has also the following conse- 
quence: 

Proposition 4.2. Let (ao, . . . , an-2, On.-i) G B"^^ xT be independent ran- 
dom variables with rotationally invariant distribution. Then 

GW(ao)e(i)(ai) . . . G("-i)(a„_i) '^= H(°)(ao)H«(«i) ■ • • H("-i)(a„_i). 

Proof. We give two proofs of this result. The first one is a consequence of 
Theorem 12.81 from which we know that 



e(°)(ao)e(i)(ai) . . . e("-i)(a„_i) = hW(7o)h(i)(7i) . . . E^"-'\^n-i) , 

and the remark at the beginning of Section [3l 

For the second proof, we proceed by induction on n. For n = 1 the result 
is obvious. Suppose the result holds at rank n — 1 : thanks to the recurrence 
hypothesis, 

HW(ao)H«(ai) . . . H(-i)K_i) 'i=" HW(ao)eW(ai) . . . e(--i)K_i). 
Let e^^" = . ~^" . An elementary calculation gjives 

HW(ao)e(i)(ai) . . . e("-2)K_2)e("-i)(a„_i) 

As the Ofc's are independent with law invariant by rotation, 

(qo, e""^"ai, . . . , e~"^°an-2,e"^°an-i) = (ao, ai, . . . , an-2, a^-i), 
which completes the proof. D 

Now that we have a matrix model for the circular Jacobi ensemble, we 
can study the characteristic polynomial for such matrices; the key formula 
will be (IXTT|) . 

Proposition 4.3. Let U be a unitary matrix of size n and let Z„ = det(Id — 
U) be its characteristic polynomial evaluated at 1. Then, in the circular 
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Jacobi ensemble, Zn can be written as a product of n independent complex 
random variables: 

n-l 

Zn='[\{l -Ik), 
k=0 
where the laws of the 7^ 's are given in Theorem \3.^ Consequently for any 
s,t € C, with 9^e(t) > — 2, the Mellin- Fourier transform of Zn is : 

lio r(/3'^ + 1 + 6+ ~5)T{^'k + 1 + 5 + ^)r(/3'A; + 1 + 5+'-^)' ^ ' ' 

Proof. The first part is an easy consequence of (|2.1ip and Tlieorem 13.21 To 
prove tfie second part, we note that if X^ = (1 — 7fc), then |Xfc|*e'''^'^§ *= = 
(1 — 7fc)"(l — 7fc)^) where a = {t + s)/2 and h = {t — s)/2. Consequently, by 
independence of the 7a:'s, we obtain: 

n — 1 

nZnfe^'^'^^-] = n IE[(1 - lkT{l - lk)\ 
fc=0 

and formula (j4.2p then easily follows from Lemma |3.1[ D 

5. Limiting spectral measure and large deviations 

In (jl.7p we defined the spectral measure which is a central tool for the 
study of our circular ensembles. Let us re- write this measure on T as 

n 

fc=i 

where we put a superscript *•"' to stress on the dependency on n. Besides, 
in classical Random Matrix Theory, many authors are mainly interested in 
the empirical spectral distribution (ESD) defined by 



/^e- = ^E^4"). (5.2) 

'^ fc=l ^ ' 

(n) 

We are concerned with their asymptotics under C J^ J when n — ?■ 00 with 

6 = 6{n) = /3 nd 

where IHe d > (and as usual /?' = /3/2). In this framework the variables 

{9i , . . . , On ) and (vrj"" , . . . , 7r„ ) are distributed as in Theorem 13.31 It is 

well-known that in the CUE (i.e. for d = and /?' = !), the sequence (/iggd) 
converges weakly in probability to the uniform measure on T and that it 
satisfies the Large Deviation Principle at scale n^ ([T7| chap. 5). 



In this section we prove that in the cJUE, both sequences (/Xsp ) and (/ig' 



(")^ ^^A ^,(")^ 



sd/ 



converge weakly in probability to the same limit measure supported by an 
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arc of the unit circle, and that (/igsd) satisfies the Large Deviation Principle 
at scale nP. In the next two subsections, we first recall some definitions 
(which may be skipped by probabilist readers) and then state the main 
results which are Theorems 15. H 15.41 and 15.51 



5.1. Weak convergence of measures. Let 7Wi(T) be the space of prob- 
ability measures on T. The weak topology on 7Wi(T) is defined from the 
duality with the space of continuous functions on T; that is, u^ weakly con- 
verges to i^ means J fdvn -^ J fdv for every continuous function /. Since 
T is compact, it is equivalent to the convergence of moments. If we define 
the distribution function oi v ^ Mi{T) as F^{t) = iy{{e'^;9 E [0, i)}), then 
i^n ^ i^ weakly if and only if 

for all t € [0, 27r) at which F,^{t) is continuous. The Levy distancqj d-L is 
compatible with this topology and makes A4i(T) a compact metric space. 
It should be noticed that for all pair of elements of A^i(T) 

dL{^^,u)< sup \F^{t)-F^{t)\. (5.3) 

t 

In the following we consider random probability measures (f„) and the 
weak convergence in probability of (f^) to a deterministic u £ Ali(T). That 
means that for every n > 1 we have a probability space {i}n,Fn,Qn), a 
A^i(T)-valued random variable f„, and that d{un,i') — ?■ in probability, 
where d is any distance compatible with the weak topology (e.g. the Levy 
distance d^). In other words 

limf„ = 1/ (in probability) <^ lim (i(z/„, z^) = (in probability) 

n n 

^ Ve > lim(5„(d(i/„,i^) > e) = 0. 

n 

(5.4) 

5.1.1. The spectral measure. The following theorem states the convergence 
of the sequence (/^sp ) to an explicit limit, which we identify now, with the 
help of some more notation. We assume that 9^ed > 0. Let 

"d = --— ^ (5-5) 

1 + d 

and let 6*^ € [0, vr) and .^d £ [~l^d|) |^d|] be such that 



. C'd 

sm — 
2 



d 



1 + d 



1 + d 



whose precise definition is not needed here (see [10] Theorem D8) 
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If 9^e d > 0, d / 0, let Wd be defined by 



Wd{e) 



|l + ad| sin(6'/2) 
otherwise. 



In tliis case let 
and for d = let 



(5.6) 



riff 

d^,d{C)■■=wd{e)^ iC = e''), (5.7) 

d/^o(C):=^ (C = e'^), 



the Haar measure on T. 

Theorem 5.1. Assume D\e d > 0. As n — t- oo, 

lim n^^' = /id (in probability) . 

To prove this convergence we use the parameterization of measures by 
their modified Verblunsky coefficients and the following two lemmas whose 
proofs are postponed until the end of the section. 

Lemma 5.2. For every fixed k >0, as n —^ cxo, 

lim7^" = = (in probability) , (5-8) 

and consequently 

lima^ = Odc"'^ '^'^ (in probability) . (5-9) 

Lemma 5.3. The sequence of Verblunsky coefficients of ^^ is precisely 

Proof of Theorem 5.1. For v G Mi{T), let mk{v) = f e'^'^du{9) be its k-th. 
moment and let afc(z^) be its A:-th Verblunsky - in short V- coefficient. Mo- 
ments are related to V-coefficients in a continuous way : for every j > 1, 
m,j{v) is a continuous function of (qq, . . . , a(jAAf)-i) where A^ is the cardi- 
nal of the support of v. We know for Lemma 15.21 that for every fixed A:, 
limaj(//sp ) = «j(^d) for j < A:, in probability. It entails the convergence 

of mfc(/isp ) to mfc(//d); in probability, hence the weak convergence of the 
measures. D 

Remark. The above method can easily be adapted to show the trivial asymp- 
totics of two different scaling regimes : 

• if 5{n) = o{n), the limit is the uniform measure on T, 

• if 5{n) is real and n = o{5{n)), the limit is the Dirac measure at —1. 
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Proof of Lemma \5.2[ For 7^" we use the Mellin transform 

E^fl- (n)^ A _ r(/3^(n -k-l) + 6 + 6_ + s + l)r(/3^(n - fc - 1^ + ^ + l) 
V ^^ ' J T{P'{n-k-l) + 5 + 6 + l)T{l3'{n-k-l)+5 + s + l)' 

(this comes immediately from (13. 3p ). 
Since for fixed z £ C 

n.-s>oo r(n)n^ 
we get, for fixed s (and A;) 

V in //I {n)^s^ /1 + d + d 



1 + d 
which implies that 

lim (1 — 7^," ) = = — in probability . 

n^oo " 1 + d 



which is equivalent to (I5.8p . The statement (]5.9p is a direct consequence of 
(USD and (I22D. D 



Rem,ark. The convergences in probability in the above lemma actually hold 
in LP, for all p > because all variables are bounded by 1. 



Proof of Lemm.a 15.31 If d is real, the Verblunsky coefficients are all equal. 
In fact, measures satisfying this property are known to have an absolute 
continuous part supported by an arc of the unit circle, and a possible ad- 
ditional Dirac mass ([35j p. 87). More precisely, let a € D and let v be the 
measure on T such that afc(i^) = a for every k > 0. If 6{a) and ^{a) are 
defined by 

e(a) = 2arcsin \a\ , e'«(") = ^-^ , 

1 + a 

then 

rif) 

du{C) = w{9)- (C = e'^), 

where 



w{e) 



sm\e/2)-sm\e{a)/2) 

l + a|sin((g + g(a))/2) if ^ ^ W«)> 2. " ^(«) , (5.10) 
otherwise 



as soon as 



l« + ^l<^ (5.11) 

(otherwise, there is an additional Dirac mass). The orthogonal polynomials 
with respect to this measure are known as the " Geronimus polynomials" . 
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When a = a^, we set 6{a) = 6d, C(q^) = '^d and v = Vd, w = Wd- We see 
that _ 

1 1-d 

""^+2 = ^(171)' 
so that the condition (jS.lip is fulfilled if and only if IHe d > 0, which we 
assumed. 

Moreover, it is known (see [3^ p. 960) that if {ak)k>o is the sequence of 
Verblunsky coefficients of a measure //, then the coefficients {e~^^'^~^^'^'^ak)k>o 
are associated with /j, rotated by S,d- Consequently, 

dfidiO = dudiCe-'^'^) , 

which is precisely (j5.7p . D 

5.1.2. The BSD. In matrix models, the convergence of the ESD is often 
tackled directly via the convergence of moments or the convergence of the 
Cauchy transform. Here, we follow a different way: we use Theorem 15.11 and 
prove that the two sequences (/Xsp )■« and {^Jiesd)n are "contiguous". 

Theorem 5.4. Assume 9^e d > 0. As n — )■ oo, 

lim/igg^ = /Xd (in probability) , 

n 

where /id is given in (5.1). 

Proof. Prom Theorem 5.1, we know, with the notation from the beginning 
of Section 15.11 that 

lim (ii(/i^",'', /id) = in probability . 
By the triangle inequality 

rfL(/ii:L/^d) < dL(/ig\/^d) + dL(/i(;;),/.ii) 

so, it is enough to prove that ^^(/isp , /^esd) converges to in probability. 
Thanks to (|5.3p . it is enough to prove 

lim sup |F („) (t) — F („) {t)\ =0 in probability . (5-12) 

n— )-00 J Msp Mesd 

In fact, 

sup |F M (t) - F („) (t)| = max |4"^ - -| , (5.13) 

where S^ = '^j^i vr:"' . Using the union bound and the Markov inequality, 
we may write 

k 



I max 



q(«) 




n 



> e 



^)1 • f^'"' 
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and then use explicit distributions. Indeed, we showed in Theorem 13.31 that 
the vector (vrj" , . . . ,7rn ) fohows the Dir„(/3') distribution. It entails that 
for A: = 1, . . . , n — 1, the variable S^ is Beta(/3'A;, /3'(n — k)) distributed. 

Recall that the Mellin transform of a beta variable Beta(a, b) with positive 
parameters a and b is 

E(Beta(„.^.)= ^i°;;)^'° + ^i . 

r(a)r(a + b + s) 

for s > —a. As a consequence, EBeta(a, b) = a/{a+b) and a straightforward 
calculation gives 

3a6(2a2 + 26^ - 2a6 + a'^b + afc^) 
(a + 6)4(a + 6 + l)(a + 5 + 2)(a + 6 + 3) 

o- "^ 



E (Beta(a,6) - EBeta(a, 6))) 



When a = 13' k and b = f3'{n — k), this shows that the k-th term in the sum 
of (l5llDisOf^^),sothat 

>e]=0\\ ^^^^^]=o(^ 



I max 



'^ n 




Thanks to (j5.13p . this yields (j5.12p and completes the proof. D 

5.2. Large deviations for the ESD. It turns out that the convergence in 
Theorem 15.41 is exponentially fast, and to make this statement more precise 
we need to recall the definition of some notions of large deviations, for 
which the reference is the book of Dembo and Zeitouni ([lOj). We say that 
a sequence (P„) of probability measures on a measurable Hausdorff space 
{X,B[X)) satisfies the LDP at scale Un (with Un —> oo), if there exists a 
lower semicontinous function I : X ^- [0,oo] such that 

limsup — logP„(F) < -inf{/(x);x € F} (5.15) 

Un 

liminf — logP„(G) > -inf{/(x);x € G} , (5.16) 

Un 

for every closed set F C X and every open set G C X. The rate function 
/ is called good if its level sets are compact. More generally, a sequence of 
^-valued random variables is said to satisfy the LDP if their distributions 
satisfy the LDP. Prom now on, we work with X = A^i(T) whose compacity 
makes our task simpler. In fact, according to |10j Theorem 4.1.11, in this 
case the LDP is equivalent to the following property: for every fi € 7Wi(T) 

— I{fJ.) = limlimsup — log-P„(i?(^,e)) (5-17) 

si-O n Un 

= limliminf — logP„(S(^,e)) (5.18) 

e;0 n Un 
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where B{fi,e) is the open ball of radius e, centered at n 

B{^i, e) = {u : diifJ., u) < e} . 
It is easy to see that the latter is equivalent to the pair of inequalities 

limlimsup — logP„(S(;U,e)) < -I{fi) (5.19) 

limlimsup — logP„(S(;U,e)) > -I{n) . (5.20) 

e4,0 n Un 

For the sake of completeness, let us explain shortly why (|5.19p and (j5.20p 
lead to (J5.15P and (j5.16p . respectively. From the one hand, every closed 
(hence compact) set in Ali(T) may be recovered by a finite number of balls 
and applying ()5.19p for all these balls leads to ()5.15p . On the other hand, 
every open set contains open balls and it remains to optimize with repect 
to the centers of the balls. 

Our large deviations result follows the way initiated by the pioneer paper 
of Ben Arous and Guionnet ([5]) and continued by Hiai and Petz (p^.|18]). 

We work with the set A^i(T) of probability measures on the unit circle. 
For fi G Ali(T), the Voiculescu entropy is 

S(^) = J J log IC - C'\dKC)dKC') . (5.21) 

With a different sign it is the logarithmic energy of /x (see Chapter 5.3 in 
[It]). We also define the potential 

' -2{me d) log (2 sin |) - (Jm d){e - tt) if C = e'^ 9 € (0, 27r) 
oo if C = 1 and 9^e d > 

-prndlTT if C = 1 and 9^ed = 

(5.22) 



QdiC) 



It should be noticed that Qj is a lower semi continuous function (l.s.c). 
The main result of this section is an extension of the theorem of Hiai and 
Petz (Theorem 5.4.10 in [17J), which corresponds to the case d = + a 
continuous potential. 

Theorem 5.5. Let d e C with 9^e d > 0. For n e N and (Ci, • • • , Cn) G T" 
let 

n 

MCl, • • • ,Cn) := 11(1 - ar'-ii-ar''^ n 10 - ai'^' 

fc=l j<k 

(n) 

and let P^ be the distribution on T" having the density 

/i(Ci,--- ,Cn), (5.23) 



where -Zd(n) is the normalization constant. 
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(1) We have 

lim 4^1ogZd(n)=S(d) (5.24) 

where 

B{d) = I [x\ogx+{x + 29^e d) log(x + 29^e d)] dx 
Jo 

— I \{x + d) log(x + d) + (x + d) log(x + d)] dx . 
Jo 

(2) The sequence of distributions of 

,M _ ^1 + • • • + (^Cn 

Mesd ~ 



n 

under F^ satisfies the LDP in Ali(T) at scale /3'n^ with (good) rate 
function 

hifi) = -S(^) + / QdiOdKC) + B{d) . (5.25) 

(3) The rate function vanishes only at fi = fid- 

Proof. (1) An exact expression of 2{n) is obtained using the following 
lemma, whose proof is postponed to the end of this subsection. 

Lemma 5.6. // we define 

„ n 

ZsM = 17(1 - e^'^ni - e-''^Y n \e''^ - e^'^fdOi . . . dO^, 

then we have 

^ , ,_ Ti/3'n + l) ^^ T{f3'j + l)T{f3'j + l + s + t) 

^^'*("^ - (r(/3' + 1))" y r(/3'j + 1 + s)TW'j + 1 + 1) • ^'-'^^ 

We have Zd{n) = ^d/3'nd/3'n(^) ^^^ then, taking for log the principal 
value of the logarithm, 

logZd(n) = logr(/3'n + l)-nlogr(/3' + l) 

n— 1 

+ J2 N r(/5'i + 1) + log ^(^'J + 1 + 29^e d/3'n)] 

j=0 

n— 1 

- Y. N r(/5'i + 1 + d/?'") + log r(/3'i + 1 + d/3'n)] 

j=0 

+ 2iA;„7r (fc„ € Z , \k\ < 5n) . 



30 



P. BOURGADE. A. NIKEGHBALI. AND A. ROUAULT 



From the Binet formula (Abramowitz and Stegun pLj or Erdelyi et al. [2] 
p. 21), we have for IHex > 

I 1 r°° 

logr(x) = (x- -)logx-a; + -log(27r)+ / f{s)e-'''ds. (5.27) 

where the function / is defined by 

1_2^ 1 



1 


1 1 1 





\ 7 


2 


s e** — 1 



k=l 



s2 + 47r2fc2 ' 



f{s)- 
and satisfies for every s > 



< fis) < /(O) = 1/12 , < Isfis) + 2 ) < 1 • 
Using ()5.27p . a straightforward study of Riemann sums gives 



1™ l^m ^°S ^d(n) = S(d) 



where 



5(d) = / [2;logx + (x + 2$Hed)Iog(x + 2fRed)] dx 

-/ [(x + d)log(x + d) + (x + d)log(x + d)] dx. (5.28) 
Jo 

(2) The proof is based on the explicit form of the joint eigenvalue density. 
We follow the lines of [S], [IS], [IS] and [U]. We skip the index esd for 
notational convenience. Our goal is the proof of the two inequalities (cf. 
(f5T9]l and dOOD ). which hold for every fj. € Mi{T): 

limlimsup— ^logPj'^(^(") € B{n,e)) < -/d(/u) (5.29) 



e^O 



»(") (-,,(") 



limliminf — ^logPr(^W ^ B{fi,e)) > -/d(/u) 



(5.30) 



We rest on the LDP known in the case of d = 0, and use the general 
method which consists in estimating the Radon-Nikodym derivative ruling 
the change of probability. The key formula is 



Pi'^)(/.(")EB(/.,e)) = |44<^ 



'^^i(")eB{^l,e)'^ 



-n'^13' jQddfi^"') 



(5.31) 



The upper hound (15.29)) . Let us first assume 9^ed > 0, so that QdiC) ~^ cxd 
as C ^ 1- For i? > we consider the cutoff Q^ = min(Qd5 R)- Since Q^ is 
continuous, the mapping v € A^i(T) H' f Q^ du is continuous, so 



inf [Q^du> [ Q^dfi-ri{e,R) 



(5.32) 
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with lim^ ri(e, R) = 0. Since Qd > Q we get 



LlogE^ 



/3'n 



l^(")e-B(At,e)^ 



-n^P' fQadti<-"^ 



< 



1 



- Q^dfi + ri{e,R) (5.33) 

Thanks to the LDP known for d = 0, we can take Umsup in n and then hmit 
in e and obtain 

lim hm sup log Ei" 



P'n^ 



1 , , p-n^P'jQidf,(") 



<S(/i)- / Q^d^ 



^(")e_B(/i,£)^ 
and thanks to (fOT]) and (fOD 

limnmsup-^logP!,"^(/i(") eB{fi,£))<J:{fi)- ( Q^ dix- B{d) . 

By the monotone convergence theorem hmi^^oo / Q dfj, = J Qddfi, which 
proves (fOOD . 

If IHe d = 0, then Qd is lower semi continuous and bounded, so that the 
mapping v ^ J Qddv is lower semi continuous and (|5.32p still holds with Qd 
instead of Q^ and some r{e) instead of ri{e,R). The rest of the argument 
is the same as above. 

The lower bound i5.30\) . If /d(/^) = oo, the bound is trivial, so that we 
can assume that // has no atom. To overcome the problem of the singularity 
at 1, we use a classical approximation of /i by a probability vanishing in a 
neighborhood of 1, i.e. 

The benefit is that the mapping i/ i— )• / Qddv is continuous in a neighborhood 
of fiM- Choosing M large enough, we ensure that dL{f^M, IJ') < s/2 and then, 
thanks to the triangle inequality 

B{ti,e) D B{jiM,e/2), 

which leads to 

Pi")(M(") G B{fi,e)) > Pi")(/.(") e Bi^,M,e/2)) . 

We come back to ()5.3ip with fiM and e/2. For u G B{fiM,£/2), we have 
j Qddv < j Qd dfiM + r2{e, M)) where lim^ r2(e, M) = 0. We get 



yo^<' 



I3'n 



l^(")eS(M,£)^ 



-n'^P' fQad,i<-"'> 



> 



1 



f) (/.(") Gi?(MM, e/2)) 



- / QddfJ-M -r2{e,M) . 



Again, we take liminf in n and lim in e, and use successively the LDP for 
d = 0, ([OT]l and ([Oijl to obtain: 

limlimsup-i^logPi"V(") G B{fi,e)) > S(/iM) - [ QddfiM - B{d) . 
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Now, since I{iJ:) < oo and since Qd and S are bounded below, the monotone 
convergence theorem yields 

lim T.{fiM) = Y,{fi) , lim QddfiM = Qddfi. 
Af — >oo _A/-j-oo J J 

This ends the proof of (j5.30p . 

(3) The uniqueness of the minimizer is a direct consequence of the strict 
convexity of I which comes from the strict concavity of S. 

We do not give a self contained proof of the identity of the minimizer, 
but rather use a probabilistic argument. On the one hand in Theorem [5 



(n) 

we proved that /ig^^ converges weakly in probability to /id, and on the other 
hand the LDP combined with the uniqueness of the minimizer imply that 

(n) 

Mesd conv 
proof. D 



(n) 

Mesd converges weakly in probability to this minimizer. This completes the 



Proof of Lemma \5.6[ We have 



2o,o{n) 



E(det(Id-C/)Met(Id -[/)*) 



\t 



(n) 

where the mean is taken under the C Jq i distribution. Under this distribu- 
tion, det(Id — U) has the same law as the product of independent variables 
1 — Ok, where Ofc is z/^(„_fc_i)+i distributed. We get 

E(det(Id - C/)Met(Id - C7)*) = [J E(l - ajYil - aj)Y 

j=0 

From (jS.Sp we get 

Zo,o(n) y r{(3'j + 1 + s)ril3'j + 1 + t) 



Besides, Lemma 4.4 in [21] gives 

T{P'n + l) 
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